Phenomenological insight into JLab proton polarization data puzzle by deuteron 

impulse approximation 
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The non-relativistic impulse approximation of deuteron electromagnetic form factors is used to 
investigate the space-like region behavior of the proton electric form factor in regard of the two 
contradictory experimental results extracted either from Rosenbluth separation method or from 
recoil proton JLab polarization data. 
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I. INTRODUCTION 



A study of the deuteron electromagnetic (EM) struc- 
ture is very interesting as the deuteron is the most sim- 
ple bound state system of nucleons, which provides an 
excellent opportunity to study nucleon-nucleon (NN) in- 
teraction as well as its dependence on electromagnetic 
structure of underlying nucleons. Up to now several (phe- 
nomenological) models and fits of the deuteron electro- 
magnetic structure have been developed. The purely phe- 
nomenological fits [l[ give small % 2 but they come with- 
out any physical background. The vector meson domi- 
nance (VMD) models (and their generalization by incor- 
porating correct deuteron form factor (FF) analytic prop- 
erties) Hii ii describe deuteron EM structure through 
exchange of isoscalar vector mesons u, <p (and their exci- 
tations) and they don't assume NN interaction effects ex- 
plicitly. Another class of models, non-relativistic (NIA) 
and relativistic (RIA) impulse approximations, reviewed 
in assume NN interaction based on deuteron wave 
functions and describe the deuteron EM FFs through 
the nucleon EM FFs. Just such models seem to be able 
to bring some light into existence of two contradicting 
space-like behaviors of the electric nucleon FF, obtained 
in elastic scattering of unpolarized electrons on unpolar- 
ized protons by Rosenbluth method and in polarization 
transfer process e~p — > e~p measuring transversal and 
longitudinal components of the recoil proton's polariza- 
tion simultaneously. 



II. TWO CONTRADICTING PROTON 
ELECTRIC FORM FACTOR BEHAVIORS IN 
SPACE-LIKE REGION 



The electron-proton elastic scattering was used to be 
the most common way of the proton electromagnetic 
structure study from the half of 50 's of the last century 



|6| and abundant data (more than 400 data points) on 
proton electric GE P (t) and magnetic GjWp(t) FFs in the 
space-like region — Q 2 = q 2 = t < appeared (for ref- 
erences see paper |7[). They have been obtained from 
the measured differential cross-section of the elastic scat- 
tering of unpolarized electrons on unpolarized protons in 
laboratory system 



da lab (e~p -> e~p) _ a 2 cos 2 (6>/2) 
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.[A(t)+B(t) tan 2 (8/2)], (1) 
a = 1/137, E-the incident electron energy 
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by the Rosenbluth technique. Their ratio 

fJ> P G Ep(t) I 'G Mp{t) in error bars roughly equals one, 
showing the electric and magnetic distributions in the 
proton to be equal. 

Recently in JLab measuring transverse 



P t = - (-2) Vt{1+t)G Mp Ge p tan(0/2) (3) 
Jo 

and longitudinal 



Pi 



h(E + E') 
I m p 



v / r(l + r)Gi fp tan 2 (0/2) 
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components of the recoil proton's polarization (as sug- 
gested in Refs. |ll|) in the electron scattering plane of the 
polarization transfer process e~p — > e~p (h is the electron 
beam helicity, Iq is the unpolarized cross-section exclud- 
ing OMott and t = —t/Am 2 ) simultaneously, the very pre- 
cise and surprising data on the ratio /i p G Ep(t) / 'G Mp(t) 
have been obtained, showing the electric and magnetic 
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distributions in the proton to be different, contrary to 
what was followed from Roscnbluth data. 

The latter contradiction is now well known as the JLab 
proton polarization data puzzle and a natural question is 
arisen "Which data are correct"? 

A lot of effort has been devoted to the search of a 
definite solution [l2j]-[l!|. The explanation is likely to be 
related to radiative corrections, which may reach 40% on 
the unpolarized cross-section, but are negligible on the 
polarization ratio. 

In Ref. [l2j the structure functions method was ap- 
plied for an evaluation of the radiative corrections in 
transferred polarization experiment and their size was 
of order of magnitude too small to bring the polarization 
data in agreement with the Rosenbluth ones. 

In the papers [HI], [HI authors came to a conclusions 
that the presence of the two photon contributions could 
solve the problem, however, the kinematical properties 
related to fast decreasing form factors and the possi- 
ble presence of inelastic contributions in the intermediate 
state were not investigated in detail. 

The structure functions method has been also applied 
in Ref. [l5j] , where it was shown that the corrections can 
become very large, if one takes into account the initial 
state photon emission. However, the corresponding kine- 
matical region is usually rejected (by appropriate selec- 
tion of the scattered electron energy) in the experimental 
analysis. 

In the paper [l(| a model dependent calculation based 
on GPD showed an agreement more qualitative than 
quantitative with one set of data [l7l ]. 

In Ref. [18] it was noted that the reason of the dis- 
crepancy lies in the slope of the reduced cross-section as 
a function of the virtual photon polarization. 

In the paper [191 ] higher order radiative corrections for 
initial state emission, in frame of the structure function 
approach, has been applied to polarized and unpolarized 
cross-section. It was shown that they have the right size 
and behaviour to bring data in agreement. The 27 pho- 
ton box is shown to be irrelevant to this problem in the 
approximation of rapidly decreasing FFs and an equally 
shared momentum between the two photons. 

In this work we take another approach, a phenomeno- 
logical point of view to investigate this problem. We show 
that the proton electric FF GEp(t), and in no case the 
magnetic FF, is responsible for the discrepancy in the 
above-mentioned sets of results. 

For the latter there are indications already in the form 
of the differential cross-section ([lj, where the proton 
magnetic FF is multiplied by —t/(4:m 2 ) factor, i.e. as —t 
increases, the measured cross-section (1) becomes domi- 
nant by G 2 Mp (t) part contribution, making the extraction 
of G Ep (t) more and more difficult. As a result, one can 
have confidence only in the proton magnetic FF data ob- 
tained by the Roscnbluth technique in higher values of 
momentum transfer squared. 

A definite conclusion for the above-mentioned asser- 
tion follows from the global analysis [2(|, [2l| of the nu- 



cleon FF data in the framework of the ten-resonance uni- 
tary and analytic model of the nucleon EM structure 0] , 
which provides a test of: 

• consistency of the JLab proton polarization data 
with all other proton and neutron EM FF data 

• consistency of these data with the powerful tool of 
physics - the analyticity 

• their consistency with the asymptotic behaviour as 
predicted (up to logarithmic corrections) by QCD. 

In the analysis all Rosenbluth data for proton elec- 
tric FF in space-like region were substituted for the ra- 
tio H P GEp(t)/GMp{t) data obtained in JLab proton po- 
larization experiments. As a result, the parameters in 
comparison with those in [3] were changed very little, a 
perfect description of the new JLab data was achieved 
and almost nothing has been changed in the descrip- 
tion of GM P (t), GEn(t), GMn(t) in both the space- like 
and time-like regions and |Gsp(t)| in the time-like re- 
gion. However, the new JLab data strongly require an 
existence of the zero (see the full line in Fig. la), e.i. a 
diffraction minimum in the space- like region of |G£ P (i)| 
around t = -Q 2 = -13 GeV 2 . 

Consequently, there are two contradicting behaviors 
(see Fig. la) of the proton electric FF GE P (t) in the space- 
like region. In the following sections they are tested in 
the framework of the NIA by existing deuteron electro- 
magnetic structure data. 

III. NON-RELATIVISTIC IMPULSE 
APPROXIMATION FOR DEUTERON EM 
STRUCTURE 



In the one-photon exchange approximation the 
deuteron EM structure can be described by three scalar 
functions - EM FFs [22j and the matrix element of the 
deuteron EM current can be written, in the most general 
form, as 

- (d' \J% M \d) = G 1 (q 2 ){t 4)dp + G 2 (q 2 ) lUe-q)- 

- e;(C q )]~G 3 ( q 2 ) ^ q)iC *- q) d,, (5) 

where £, £' are polarization vectors of incoming and out- 
going deuterons d, d' of four- momenta p^ , p'^ obeying the 
relations 

Z'.p' = ; £.p = ; = -1 ; e = -1 ; d il =p' ll +p ll ; 

However, from practical point of view another linear com- 
binations of deuteron EM FFs are used [5| 

G c (q 2 ) = G^ + ^Ggiq 2 ) 

G M (q 2 ) = G 2 (q 2 ) (6) 
G Q (q 2 ) = G^q 2 ) - G 2 (q 2 ) + (1 + rj)G 3 (q 2 ) 
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Fig. 1: Results of the analysis of Rosenbluth (dashed lines) and JLab (full lines) data by nucleon EM structure U&A model. 



where r\ = -^r, Gc(q 2 ) is deuteron charge form factor, 

Gif(g 2 ) is deuteron magnetic form factor and GQ(q 2 ) is 
deuteron quadrupole form factor. 

The calculation of deuteron EM FFs within impulse 
approximation requires a knowledge of the deuteron wave 
function and nucleon EM FFs. As deuteron can be 
found in S- (» 96%) and D-state (« 4%), then NN non- 
relativistic full wave function of the deuteron can be writ- 
ten in terms of two scalar wave functions 

(/, 1, 771 — m s , m s |l, m) = 

= ^Iyo^xI? (7) 
r 

+ ^2 Y 2< m-m,xlT e (2,l,m-m s ,m s |l,m), 

where (1, 1, m — m s , m s \l, m) are Clebsh-Gordan coeffi- 
cients, Yi_ mi are spherical harmonics normalized to unity 
on the unit sphere and zq = u, z-i — w are reduced S— 
and D-state wave functions, respectively. 
The normalization condition 

/ d 3 r^ abm ,^ abm = S m , m 



implies normalization 

POO 

/ dr [u 2 (r) + w 2 (r)] =1, (8) 
Jo 

which could be understood as probability of finding 
deuteron in S— or D-state. The D-state probability 

/>oo 

Pd = drw 2 (r) 
Jo 

is an interesting measurement of the strength of the ten- 
sor component of the NN force. 

The best non-relativistic wave functions are calculated 
from the Schrodinger equation using a potential adjusted 
to fit NN scattering data for lab energies from to 350 
MeV. In this paper we will use one of the most common 
potentials called Paris potential [231 ] . which was among 
the first potentials to be determined from such realistic 
fit. The S— and D— state wave functions determined 
from this model are presented in Fig. O 

The deuteron is an isoscalar target, therefore within 
non-relativistic impulse approximation, its FFs depend 
only on the isoscalar nucleon form factors G S EN and G S MN 

Gen = Gep + Geu 

Gmn = Gmp + Gmu (9) 
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Fig. 2: The S— and D— state wave functions (u and w, re- 
spectively) behaviors for Paris potential. 



in the following way 



Gc — G s EN Dc 
m d 



Gm 



2m Tl 



[G S MN D M + G S EN D E ] 



(10) 



Gq — G s en Dq 



where the body form factors Dc, Dm, De and Dq are 
functions of the momentum transfer squared t. The non- 
relativistic formulas for the body form factors D involve 
overlaps of the wave functions u(r),w(r), weighted by 
spherical Bessel functions 



Dc(q 2 ) 
D M {q 2 ) 

D E (q 2 ) 
D Q (q 2 ) 



dr [u 2 {r) + w 2 (r)] j (n) 

dr[2u 2 {r)~w 2 {r)] Jo ( K ) 
o 

V2u(r)w(r) + w 2 (r) J2(k) 

drw 2 (r) [j'o(k) + J2(k)] 
w(r) 
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(11) 



V2r]Jo 



drw{r) 



u(r) 



J2(«) 



where k = qr/2. At q 2 — 0, the body form factors be- 
come 



D c (0) 
D M (0) 
D e (0) 
Dq(0) 



dr [u 2 (r) + w 2 (r)] = 1 

/ dr [2u 2 (r) - w 2 (r)] =2- 3P D 
Jo 

2 J drw 2 (r) = -P D (12) 
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Fig. 3: Deuteron structure function A(t). 
giving the non-relativistic predictions 

fi s N (2-3P D ] 



/'<J 



Dq(0) 



(13) 
1-5P D , 



where Qd is the quadrupole moment of the deuteron, 
fid is the magnetic moment of the deuteron and fj, s N = 
i(/i p + fi n — 1) is the isoscalar nucleon magnetic mo- 
ment. The experimental value of the deuteron magnetic 
moment [id — 1.7139, leads to probability of D-state 
Pd = 4.0%. But this is only approximate value, because 
the magnetic moment is very sensitive to relativistic cor- 
rections. 

Experimentally the EM structure of the deuteron 
is measured in the elastic scattering of electrons on 
deuterons, described by the differential cross-section ((T|) 
with 



A(t) 
B(t) 



G 2 c (t) + ^G 2 M (t) + - V 2 G 2 Q (t), 



rtl + vVG 2 M (t), 



(14) 



and applying the Rosenbluth technique. As a result 
the data on structure functions A(t), B{t) are obtained, 
which are found to be compiled in the paper 



IV. RESULTS 

In order to test the two contradicting behaviors of 
GEp{t), as shown in Fig. la, in comparison with the 
deuteron structure functions A(t), B{t) data, we use 
expressions for deuteron EM FFs (fTU|) to be expressed 
through nucleon EM FFs. First, the fits of Rosenbluth 
data were made within Unitary and Analytic model of 
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Fig. 4: Deuteron structure function B(t). 



Fig. 6: Magnetic deuteron FF Giw(t). 




t [GeV] 

Fig. 5: Charge deuteron FF Gc(t). 



nucleon EM FFs [7j, then the fits of JLab proton po- 
larization data with all other existing nucleon EM FFs 
data. From both behaviors the isoscalar nucleon FFs 
were determined, by means of which the deuteron FFs 
have been found and as a result the two differen behav- 
iors of deuteron structure functions A(t), B{t) were cal- 
culated. 

By a comparison of the latter (see Figs. 3 and 4) with 
existing deuteron structure functions data corresponding 
X 2 s are calculated as presented in Table HI 

from where one can see immediately that the behaviors 
of A(t), B(t), obtained by means of the GE P (t) with the 
zero around t — —13 GeV 2 unambiguously is preferred. 

For completeness we present (see Figs. 5, 6 and 7) 




t [GeV] 

Fig. 7: Quadrupole deuteron FF GQ(t). 



also the obtained behaviors of deuteron EM FFs Gc(t), 
Gm(^) and Gq(<) with the data [j| obtained in recent 
polarization experiments. 





x\ 


Xb 


JLab 


926 


476 


Rosenbluth 


2080 


573 



Table I: The \ 2 of deuteron structure functions A(t), B(t) for 
two different scenarios. 
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V. CONCLUSIONS 

The Rosenbluth data on the proton electric FF in 
space-like region and the new JLab proton polarization 
data on the ratio \i v G E P (t) / G Mp(t) , were both analyzed 
together with all other existing nucleon FFs data in 
space-like and time-like regions by means of the ten- 
resonance unitary and analytic nucleon EM structure 
model. As a result, two contradicting behaviors oiGspit) 
in space-like region were obtained. They are brought 
into a comparison with other independent data, the data 



on deuteron structure functions, by means of the non- 
relativistic impulse approximation of deuteron EM struc- 
ture. From the values of \ 2 shown in Table I it unambigu- 
ously follows that the GE P (t) from the JLab proton polar- 
ization data analysis with the zero around t = — 13 GeV 2 
are more consistent with the deuteron structure functions 
A{t),B{t) data than the older Rosenbluth behavior. 

The work was partly supported by Slovak Grant 
Agency for Sciences VEGA, Grant No. 2/4099/26 (C.A., 
S.D. and A.Z.D.). 
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